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In this paper we discuss the origin of the asymmetry present in D meson production and its energy 
dependence. In particular, we have apphed the meson cloud model to calculate the asymmetries in 
D~ /D'^ meson production in high energy p — p collisions and find a good agreement with recent 
LHCb data. Although small, this non-vanishing asymmetry may shed light on the role played by 
the charm meson cloud of the proton. 



PACS numbers: 12.38.-t, 12.38.Bx, 13.60.Le 

I. INTRODUCTION 

It is experimentally well known [ll-Q that there is a 
significant difference between the Xp (Feynman momen- 
tum) distributions of and D~ mesons produced in 
hadronic collisions with proton, and pion projectiles. 
It is usually quantified in terms of the asymmetry func- 
tion: 
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where N may represent the number of mesons of a spe- 
cific type or its distribution in xp, rapidity y and pr- 
The recent data of the COMPASS collaboration 01 have 
confirmed the existence of charm production asymme- 
tries also in jp collisions. Moreover, the very recent data 
from the LHCb collaboration Q showed that there is 
asymmetry in the production of and D~ mesons in 
proton-proton collisions at 7 TeV. The origin of these 
asymmetries is still an open question. It is not possible 
to understand these production asymmetries only with 
usual perturbative QCD (pQCD) or with the string frag- 
mentation model contained in PYTHIA. This has moti- 
vated the construction of alternative models [9l-[l]| which 
were able to obtain a reasonable description of the low 
energy data and make concrete predictions for higher en- 
ergy collisions. The LHCb data allow us, for the first 
time, to compare the predictions of the models with high 
energy data. Moreover, studying the energy dependence 
of production asymmetries, it may be possible to learn 
more about forward charm production, which undoubt- 
edly has a non-perturbative component [loj . In this work 
we update one of these models, the meson cloud model 
or MCM [oj], and compare its predictions with the new 
LHCb data. 

Let us now briefly review some ideas about charm pro- 
duction. In perturbative QCD the most relevant elemen- 
tary processes which are responsible for charm produc- 
tion are <? + g — > c -\- c and g + g ^ c + c. At high en- 
ergies, due to the growth of the gluon distributions, the 



latter should be dominant. In standard pQCD, after be- 
ing produced the c and c quarks fragment independently 
and hence, the resulting mesons and (also 13° and 
D'^) will have the same rapidity, pt and xp distributions. 
This is indeed true for the bulk of charm production. Dif- 
ferences between the and xp distributions appear 
at large xp, with D~ being harder. Given the valence 
quark content of the proton p(uud) and of the D~(dc), 
a natural explanation of the observed effect is that the 
c is dragged by the projectile valence d quark, forming 
the, somewhat harder, dc bound state. This process has 
been usually called recombination or coalescence and it is 
illustrated in Fig. [T] This is a non-perturbative process 
and recombination models have been first proposed long 
time ago p^ - IT4| and then used more recently jH, [11] to 
study the accumulated experimental data and to make 
predictions for the RHIC collisions. Unfortunately, the 
current RHIC experimental set up did not allow for a 
precise determination of production asymmetries. How- 
ever, a more detailed analysis of the heavy quark sector 
is expected to be possible in the upgraded RHIC facility 

- RHIC H 

An alternative way to implement the idea of recom- 
bination is to use a purely hadronic picture of charm 
production, in which, instead of producing charm pairs 
and then recombining them with the valence quarks, we 
assume that the incoming proton fluctuates into a virtual 
charm meson - charm baryon pair, which may be liber- 
ated during the interaction with the target. This kind 
of fluctuation is unavoidable in any field theoretical de- 
scription of hadrons and, in fact, it was shown [31 to be 
quite relevant to the understanding of hadron structure. 
It has been also successfully applied to parti cle p roduc- 
tion in high energy soft hadron collisions [l^ l2Cll | and in 
P it has been extended to the charm sector. This mech- 
anism, in which the "meson cloud" plays a major role, 
is quite economical and can be improved systematically 
(see 2l| for light mesons). From now on we shall call 
it meson cloud model (MCM). A simple and accurate 
description of charm asymmetry production at lower en- 
ergies {^/s ~ 10 — 40 GeV) within the framework of the 
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MCM can be found in [23. 

Another popular model of forward charm production 
is the intrinsic charm model (ICM) [H IH ill. The 
existence of an intrinsic charm component in the wave 
function enhances charm meson production at large xp- 
While IC, as formulated in [Tl| is supported by sev- 
eral phenomenological analises, it is, alone, not enough 
to explain the difference between the and D~ xp dis- 
tributions. It is necessary to add a recombination mecha- 
nism (or "coalescence" ) in this kind of model to account 
for the observed asymmetries. In [Tol . [Tl| the intrinsic 
charm component of the proton is the higher Fock state 
\uudcc) and its existence is attributed to a multigluon fu- 
sion, which is not calculable in pQCD. In Q the intrin- 
sic charm component of the proton was a consequence 
of the meson-baryon fluctuations mentioned above. In 
this approach, since they "feel" the virtual bound states 
where they once were, the c and c distributions are dif- 
ferent from the beginning and when they later undergo 
independent fragmentation their difference will be trans- 
mitted to the final D mesons. Thus, in its "meson cloud" 
version, intrinsic charm may account for large x p charm 
meson production including the asymmetries. The charm 
quarks of the projectile wave function traverse the target 
and fragment independently. In this corner of the phase 
space this mechanism may be more effective than gluon 
fusion because the later generates final mesons with a, xp 
distribution peaked at zero. This possibility was explored 
in [23]. 

The appearance of the first LHCb data on asymmetries 
Q opens the exciting possibility of studying the energy 
dependence of forward charm production. Indeed, more 
than ten years from the last data on this subject, we 
have now data taken at an energy which is larger than 
the previous one by a factor of ~ 200 ! What do the 
models discussed above have to say about the energy de- 
pendence of forward charm production? What will be 
the fate of the production asymmetries ? The answer to 
this question is interesting not only to the hadron physics 
community but also to the studies of CP violation, since 
a change in the relative yield of particles and antipar- 
ticles may affect the interpretation of their decays and 
hence change the ammount of CP violation. 

Naively, we expect that perturbative processes grow 
faster and become more important than non-perturbative 
ones as the reaction energy increases. As a consequence 
the asymmetries would gradually disappear. In [23| a 
kinematical treatment of this problem arrived at the con- 
clusion that, as the energy grows, the energy deposition 
in the central region increases, baryon stopping also in- 
creases and the remnant valence quarks emerge from the 
collision with less energy. Even though the sources of 
asymmetry would still be operating, because of the in- 
creasing stopping of the valence quarks, the outgoing 
charm mesons with these quarks would be decelerated 
and their xp distribution would become invisible, buried 
under the much higher contribution from the (symmet- 
ric) central gluon fusion. In [l6j . using a recombination 




FIG. 1: D meson production in a proton - proton reaction. 
The charm pair is created by gluon fusion and then the quarks 
may fragment independently, as it happens to the c quark in 
figure, or they can "recombine" or "coalesce" with one of the 
projectile valence quarks, as it is the case of the c in the figure. 
Mesons formed by recombination are harder since they are 
"dragged" by the hard valence quark. 



approach, the authors concluded that the asymmetry re- 
mains nearly constant with energy. In the meson cloud 
approach of [22] the energy dependence of the asymmetry 
A is approximately given by: 
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where a^p is hadron-proton cross section and a^s is the 
total cc pair production cross section, which grows faster 
than the hadronic cross sections. Therefore the asymme- 
try decreases with the energy. 

In this work we use the model developed in [22I] to 
study the recent LHCb data on D'^ /D~ asymmetry and 
to check if the energy behavior of this asymmetry can 
be satisfactorily understood with this model. The paper 
is organized as follows. In the next Section, we discuss 
the asymmetry production in terms of the meson cloud 
model, presenting the main formulas and assumptions of 
the model. In Section Hill we present our results for the 
asymmetries in A and D production at SELEX (-^/s = 33 
GeV) and LHCb {y/s — 7 TeV) energies. We compare 
the results with the current data and make predictions 
for ^/s = 14 TeV. Finally, in Section HVl we summaryze 
our main conclusions. 



II. ASYMMETRY PRODUCTION IN THE 
MESON CLOUD MODEL 

A. The interaction between the cloud and the 
target 

In the MCM we assume that quantum fluctuations in 
the projectile play an important role. The proton may be 
decomposed in a series of Fock states, containing states 
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generating asymmetries. We shall refer to the first two 
processes as "indirect production" (/) and to the last one 
as "direct production" (£>). The first two are calculated 
with convolution formulas whereas the last one is given 
basically by the meson momentum distribution in the 
initial \MB) cloud state. Direct production has been 
widely used in the context of the MCM and applied to 
study n, and tt*^ production [l9| . Indirect meson 

production has been considered previously in (20j . 

Inside the baryon, in the \MB) state, the meson and 
baryon have fractional momentum j/m and j/b with distri- 
butions called fM/MBiUM) and fB/MBiys) respectively 
(we shall use for them the short notation Jm and /s). 
Of course, by momentum conservation, iim + ys = 1 and 
these distributions are related by [l^, [20] : 



FIG. 2: pp collision in which the projectile is in a \ MB) state. 
Figs, a) and b) show the "indirect" production and c) the 
"direct" D~ production. 

such as \p) = \uudcc). In the MCM we write the Fock de- 
composition in terms of the equivalent hadronic states, 
such as \p) = |I]++-D^). This expansion contains the 
"bare" terms (without cloud fluctuations), light states 
and states containing the produced charmed meson {D 
or Dg). The latter are, of course very much suppressed 
but they will be responsible for asymmetries. The "bare" 
states occur with a higher probability and are responsi- 
ble for the bulk of charm meson production at low and 
medium momentum {xp < 0.4), including, for example 
the perturbative QCD contribution. The cloud states 
are less frequent fluctuations and contribute to D pro- 
duction in the ways described below. More precisely we 
shall assume that: 

\p)^Z[ bo) + ... + |Af S) + ... + I A,^o) + IS++I?-)] (3) 

where Z is a normalization constant, \pq) is the "bare" 
proton and the "dots" denote all possible meson (M) - 
baryon (B) cloud states \MB) in the proton. The rela- 
tive normalization of these states is fixed once the cloud 
parameters are fixed. The proton is thus regarded as 
being a sum of virtual meson-baryon pairs and a proton- 
proton reaction can thus be viewed as a reaction between 
the "constituent" mesons and baryons of the projectile 
proton with the target proton (or nucleus). 

With a proton beam the possible reaction mechanisms 
for meson production at large xp and small pp (the 
soft regime) are illustrated in Fig. [21 In Fig. [5] a) the 
baryon just "flies through" , whereas the corresponding 
meson interacts inelastically producing a D meson in the 
final state. In Fig. [5] b) the meson just "files through" , 
whereas the corresponding baryon interacts inelastically 
producing a D meson in the final state. In Fig. [2] c) 
the meson in the cloud is already a D~ which escapes 
(similar considerations hold for D~ production with a 7r~ 
beam). This last mechanism is the main responsible for 



fM{y) = fB{i~v) (4) 

The "splitting function" fuiy) represents the probabil- 
ity density to find a meson with momentum fraction y 
of the total cloud state \MB). With Jm and Jb we can 
compute the differential cross section for production of 
£), Dq and Ac. In what follows we write the formulas 
for the specific case of D production but it is easy, with 
the proper replacements, to write the corresponding ex- 
pressions for Dq and Aj,. In the reaction pp ^ D~X the 
differential cross section for D production is given by: 

= $0 + + *Z5 (5) 

dxp 

where $o and <&/ refer respectively to "bare" and indi- 
rect contributions to D meson production and is the 
fractional longitudinal momentum of the outgoing me- 
son. $£i represents the direct process depicted in Fig. [5] 
c) and is given by [l^ [l^] : 

^D^—fD{xp)a''P (6) 



where fp> = and ct^p is the total pE++ cross 

section. Replacing by Ac/ Ac in Figs. 2a and 2b, 
exchanging B with M and replacing D~ by Ac in Fig. 
[5] c) we have a pictorial representation of Ac production 
in the MCM with the following expression for the direct 
process: 

<i>D = —fA{xp)a''P (7) 



where /a = Ia^/a^Do and a ^ is the total Dqp cross 
section. Analogous expressions can be written for the 
reaction tt~p DX. 



4 



B. The asymmetry 

Using ([5]), we can compute the cross sections and also 
the leading (Z?~)/nonleading(D+) asymmetry: 



made by the SELEX collaboration, the differential cross 
section for Ac production was parametrized as: 
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with riA — 2. Integrating the above expression we obtain 
the total cross section for Ac and Ac production cr^ = 
2/3 CT^, which, according to [111, can be related to the 
total c — c production cross section, CTce, in the following 
way: 



2. A 



0.1 ac 



where the last line follows from assuming $1^ = = 

$1^. This last assumption is made just for the sake of from where we finally have: 
simplicity. In reality these contributions are not equal 
and their difference is an additional source of asymmetry, 
which we assume to be less important than Since the 
"bare" states do not give origin to D' /D^ asymmetries 
(they represent mostly perturbative QCD contributions 
which rarely leave quark pairs in the large xp region) , we 



O-Q = 0.15 (Tec 



(14) 



(15) 



Using Q to obtain f\ and then inserting (O and (fT3 
into (|5]) the asymmetry A^{xp) can be written as: 



have made use of 



The denominator 



A'^ixp) 



of the above expression can be replaced by a parametriza- 
tion of the experimental data: 
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where n~ and are powers used by the different col- 
laborations to fit their data. Typically n"*" = 5 and 
n~ = 3.5, as suggested by the data analysis performed 
in [ll-Q. For the sake of simplicity we shall assume that 
n"*" = n~ = no = 5 for D mesons. Integrating the 
above expression we obtain the total cross section for 
charged charm meson production = 1/3 (T(f. As- 

suming isospin symmetry the charged and neutral (cr^° ) 
production cross sections are equal. Neglecting the con- 
tribution of other (heavier) charm states we can relate 
the D meson production cross section to the total c — c 
production cross section, CTce, in the following way: 



The behavior of ([T2]) [ Eq. ([TBI)] is controlled by the 
splitting function foixp) [fAixp)]. 



C. The splitting function 

We now write the splitting function in the Sullivan ap- 
proach [l^ ■ The fractional momentum distribution of 
a pseudoscalar meson M in the state \MB') (of a baryon 
\B)) is given by ,18, 20]: 
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where t and niM are the four momentum square and 
the mass of the meson in the cloud state and tmax = 
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(10) 



igi y/{l — y) is the maximum t, with ttib and 
TUB' respectively the B and B' masses. Following a phe- 
nomenological approach, we use for the baryon-meson- 
baryon form factor Fmbb', the exponential form: 



From the above relation we can extract the parameter 
(7(P from the experimentally measured (Tee: 



Fmbb' (t) = exp 



t- 
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(18) 



-0^ 



1.5 (Tec . (11) 

Inserting ^ and © into ^ the asymmetry becomes: 

foixp) 



where Ambb' is the form factor cut-off parameter. 

Considering the particular case where B = p, B' = 
1]++ and M = D", we insert (HT]) into to obtain the 
final expression for the asymmetry in our approach: 
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An analogous development can be performed for the 
Ac/ Ac production asymmetry, A^{xp). In the analysis 
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where 
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With the replacements gpoSa ~^ 9pDAa and a^^P — 
this same expression holds for the process pp 



(20) 



DoX. 



With the replacements gpos^ 



and (Tq 



(Tq the above expression holds for the process 
pp KcX . For the pion beam, we need also the splitting 
function of the state |7r~ > — > \D^*D~). In this state, 
the D meson momentum distribution turns out to be 
identical to (IT71) except for the bracket in the numerator 
which takes the form [— t + {{m^ — m^o. — t)/2mjja,^ ], 
and for trivial changes in the definitions, i.e., Qmbb' ~^ 
9nDD°'^ ^MBB' Fj^DD°' and Ambb' ^7tDD°*- Re- 
alizing that y — xp i'n the above equations, we can see 
that in the limit xp ^ 1, tmax —oo and the integral 
in goes to zero. In fact, it vanishes faster than the 
denominator and therefore A ^ 0. This behavior does 
not depend on the cut-off parameter but it depends on 
the choice of the form factor. For a monopole form fac- 
tor we may obtain asymmetries which grow even at very 
large xp- 

We emphasize that i.) our calculation is based on quite 
general and well established ideas, namely that hadron 
projectiles fluctuate into hadron-hadron (cloud) states 
and that these states interact with the target; ii.) our re- 
sults depend only on two parameters: A and TV. Whereas 
A affects the width and position of the maximum of the 
momentum distribution of the leading meson in the cloud 
(and consequently of the asymmetry), iV is a multiplica- 
tive factor which determines the strength of the asym- 
metry. 



Dp 



III. RESULTS AND DISCUSSION 

A. The energy dependence 

Although the recent data fS] are given in terms of the 
pseudo-rapidity rj, in order to study the energy depen- 
dence it is more convenient to use the Feynman xp vari- 
able, which may be written as: 



Xp 



2mTCosh{y) 2 mp cosh(rii) 2mTe^ 



for Tj > I. Moreover, if the transverse momentum of the 
final D meson is zero or very small, then mp ~ mp). 

In order to study the energy dependence of the asym- 
metry, we shall focus on the two energies where we have 
experimental data: ^/sl — 33 GeV and = 7000 GeV 
and construct the asymmetry ratio: 



Ra 



_ A{^) _ ^Di^)/<fp{V^) 



(21) 



Energy (GeV) 


(Tpp (mb) 


(Jcc (mb) 


33 


40 


0.04 


7000 


97 


8 


14QQQ 


110 


11 



TABLE I: Total pp cross section from ^2Q] and total cc produc- 
tion cross section from [23| as a function of the energy. The 
first two lines refer to measurements and the last one show 
model calculations described in the corresponding references. 



Using the definitions of the splitting function PT|) in © 
and then in ([5]), many energy independent factors cancel 
out and we find: 
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crPP{si) 
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(22) 



where in the last step we have used ((TT|) and assumed that 
CT^ = a^i^ P = const . aPP. Moreover, we have neglected 
the energy dependence of njj. The above ratios can be 
estimated with the recently obtained experimental data 
listed in Table HI 



With these numbers we find that increasing the energy 
the asymmetry ratio is Ra — 1/75, i.e., there is a strong 
decrease in the asymmetry. This happens because the 
responsible for the asymmetry is meson emission, a non- 
perturbative process, which has a slowly growing cross 
section. In contrast, the symmetric processes are driven 
by the perturbative partonic interactions, which have 
strongly growing cross sections. We end this subsection 
making the prediction for the order of magnitude of the 
/ D~ asymmetry in the forthcoming 14 TeV pp colli- 
sions. Using the last line of TableU setting yfs2 — 14 TeV 
and susbtituting the numbers in (1^^ we find Ra ~ 1/100 
showing the decreasing trend of the asymmetry. 



B. Predictions for the asymmetries 

Due to the lack of experimental data a direct compar- 
ison of / asymmetries in the same reaction, e.g. 
proton-proton, at different energies is difficult. However, 
we can relate and compare similar reactions. In what fol- 
lows we shall relate the two sets of data on asymmetries 
in proton-proton collisions: the low energy one taken by 
the SELEX collaboration [6] and the high energy one ob- 
tained by the LHCb collaboration Q. We will fit the 
data on Ac/A^ asymmetry. This will fix the value of 
and define a range of possible values for the cut off 
^DopAc- Using experimental information it is easy to go 
from the estimate of to the estimate of and then 
it is straightforward to calculate the associate asymmetry 
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FIG. 3: Comparison of the MCM asymmetry for Ac produc- 
tion with experimental data 0| for ^/s — 33 GeV. 
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FIG. 4: MCM prediction of the D°/D° asymmetry, Eq. ((19}. 



in Dq/Dq. In proton-proton collisions the leading charm 
mesons are those with valence quarks and hence Dq{uc) 
has the same properties as D~(dc). From this approxi- 
mate equivalence we infer the D~ /D^ asymmetry at the 
lower energy. The final step is to calculate this asymme- 
try at the LHC energy. This extrapolation can be done, 
since the model has a well defined energy dependence and 
fortunately the necessary ingredients (the cross sections) 
are already available. 



1. The coupling constants 

The coupling constant gpOA^ was estimated in several 
works with QCD sum rules [28|. Here we have chosen a 
representative number. In our analysis we will also need 
the coupling g^jj_^++^ about which, to the best of our 
knowledge, nothing is known. As a first guess we will 
then assume that 



• LHCb 
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10"^ to" 

FIG. 5: Comparison of the MCM asymmetry, Eq. IHl), with 
experiental data 8^ for D~ / . 
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FIG. 6: Comparison of the MCM asymmetry, Eq. with 
experimental data @] for D~ /D^ . Extension of Fig. [S]to the 
large xf region. 



2. The interaction cross sections of the cloud particles 

The D - proton and Ac - proton total cross sections 
in the 10 — 40 GeV range are unknown. We might know 
them better in the future from the CBM and PANDA 
experiments at FAIR. We know that the presence of the 
heavy quark makes these states much more compact than 
the proton and hence with smaller geometrical cross sec- 
tion. In the well studied case of the J/^J - proton cross 
section, a series of theoretical works gradually converged 
to a-^/'^P ~ 4 mb Moreover, it was found that, 

in contrast to the expectation of additive quark models, 
fjJ/i'P ^J/tpTT^ Inspired by these previous works we 
shall assume that: 



5pD-E++ = 9pDoA, = 5.6 



(23) 



^Dp _ ^A.p ^ 0.15 (Tff = 6mb 



(24) 
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0.15 


0.15 (Tec 


0.15 (T^f 
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TABLE II: Parameters used to calculate at ^/s = 33 GeV. 



Q 0.015 
+ 
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^ 0.01 
+ 

Q 

'r\ 0.005 
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A = 1 .95 GeV 
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FIG. 7: Prediction of the D /D'^ asymmetry for ^/s = 14 
TeV. 



3. The production cross sections of the charm particles 

The production cross section of charmed hadrons is 
measured in some cases and calculated with pQCD in 
others. In order to have an estimate of these cross sec- 
tions it is enough for our present purposes to use the 
pQCD results published in j25| . from where we have de- 
duced pT|) and p^ . All this information is summarized 
in Table [Til Starting from the definition ((20l) and using 
the numbers given in Table HIl we have: 



and also: 



N 



D 



32 TT 



32 TT 



(5.6)2 0.15 crPP 
32 TT 0.15 cr^c 



(5.6)2 0.15 crPP 
32 TT l.bCTrr 



320 



32 



(25) 



(26) 



In Fig. |3l we show the A^/Ac asymmetry in pp col- 
lisions at -y/s = 33 GeV. The experimental points are 
the results obtained by the SELEX collaboration. The 
lines are calculated with Eq. (|16l ) with the normaliza- 
tion fixed by (|25p . These data can not impose a stringent 
contraint on the model, but they can establish a range of 
acceptable values for the cut-off parameter, which values 
are shown in the figure. The outcoming numbers for A 
are those expected in this kind of meson cloud calcula- 
tion. If they had been smaller than 1 GeV or larger than 
5, this would have been an evidence against the model. 

In Fig. m we show the Dq/Dq asymmetry in pp colli- 
sions at -^/s = 33 GeV. It was calculated with Eq. ([19]) 
with the normalization factor given by ((26|). Since is 
fixed the only free parameter is the cut-off A, which, as 
it will be seen next, will be fixed so as to yield a good fit 
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FIG. 8: Prediction of the Ac/Ac asymmetry for ^ = 7 TeV. 



of the / asymmetry data from the LHCb collab- 
oration. Also in this case, the values of A used to draw 
the curves are quite reasonable. The shapes of Figs. [3] 
and 2] are correlated since they refer to the same vertex, 
where the proton splitts into a meson and a baryon 
Ac. Due to its larger mass the baryon takes most of the 
momentum and the resulting asymmetry peaks at very 
large values oixp. Complementarily, the Dq distribution 
peaks at lower values of xp^ which, nevertheless, are still 
large. The value of the cut-off A does not have to be the 
same as that used in Fig. [3] because, even though the 
coupling constant of a given vertex is always the same, 
the functional form of the form factor (as a function of 
the off-shell particle squared momentum) changes when 
the off-shell particle changes. 

Neglecting differences coming from the different 
isospin, we assume that, apart from trivial changes in 
the masses, the vertex pS++Z)^ has the same splitting 
function as the previously discussed pAc-Do vertex and 
therefore the asymmetry oi D~ production is given by 
the same expression used for the Dq. Now we try to 
reproduce the i/s = 7 TeV LHCb data with Eq. P^ . 
The only part of this expression which depends on the 
energy is the factor 7V^, which at higher energies will be 
corrected by the factor (P^ : 

^^(^/i = 33 GeV) 



iV■°(^^ = 7TeV) = 



Ra 
1 

75 



A^■°(^/i = 33GeV) (27) 



In Fig. [S]we show Eq. (HH) with 7V^(Vs = 7 TeV) and 
compare with the data, properly rewritten in terms oixp 
and with the definition ^ , which introduces a minus sign 
with respect to [s] . In spite of the large error bars we can 
see that the MCM is able to reproduce the non- vanishing 
asymmetry. The data are surprisingly sensitive to cut- 
off choices, being able to discriminate small variations in 
A. These data can not yet give a detailed information 
about the xp dependence of the asymmetry, but they 
show clearly that this asymmetry exists and also that it 
is much smaller than what we expect to find at lower 
energies for pp and than what we have already found for 
pion and S projectiles. After this close look into the data 
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points, this figure deserves a zoom out to reach the xp 
region which was scanned in previous lower energies ex- 
periments. This is shown in Fig. [6l from where we draw 
the most important conclusion of this work: the asymme- 
try definitely decreases at increasing energies, reaching at 
most 2 % at XF — 0.4. Finally, in Fig. [7] we show our 
prediction for the asymmetry to be measured at ^/s = 14 
TeV. For completeness, we present in Fig. [S]our predic- 
tions for the Ac asymmetry at ^/s = 7 TeV. 



tries in terms of a simple physical picture with few pa- 
rameters. It connects the behavior of the asymmetries at 
large xp with the charm meson momentum distribution 
within the cloud state. We can explain why we observe 
asymmetries, why they are different for different beams 
and why they decrease with increasing energies. 



IV. SUMMARY 
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